Quasiparticle mirages in the tunneling spectra of d-wave superconductors 
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We illustrate the importance of many-body effects in the Fourier transformed local density of states 
(FT-LDOS) of d-wave superconductors from a model of electrons coupled to an Einstein mode with 
energy £Iq. For bias energies significantly larger than Qo the quasiparticles have short lifetimes 
due to this coupling, and the FT-LDOS is featureless if the electron-impurity scattering is treated 
within the Born approximation. In this regime it is important to include boson exchange for the 
electron-impurity scattering which provides a 'step down' in energy for the electrons and allows for 
long lifetimes. This many-body effect produces qualitatively different results, namely the presence of 
peaks in the FT-LDOS which are mirrors of the quasiparticle interference peaks which occur at bias 
energies smaller than ~ Slo. The experimental observation of these quasiparticle mirages would be 
an important step forward in elucidating the role of many-body effects in FT-LDOS measurements. 

PACS numbers: 74.25.Jb, 74.50.+r, 74.20.-z 
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Many-body effects are known to influence the electron 
spectral function in cuprates, in particular the peak-dip- 
hump lineshape seen in the superconducting state by 
both angle resolved photoemission and tunneling 1 . The 
nature of the bosonic modes that give rise to this line- 
shape is a topic of much debate. Among the possibili- 
ties that have been discussed in the literature are certain 
optical phonons, as well as the spin resonance seen by 
inelastic neutron scattering. They all have a compara- 
ble excitation energy, Oo » 40 meV for optimal doped 
Bi2Sr2CaCu208+,5 (Bi2212), which is also near the en- 
ergy of the antinodal gap, A^. Very recently it has been 
suggested that scanning tunneling spectroscopy (STS) 
can be useful in revealing the electron-boson coupling 
in the cuprates 2 -. In STS the differential conductance 
dI/dV(r,eV) is a measure of the local density of states 
(LDOS) p{r,uj — eV) of the electrons, while the Fourier 
transform (r — > q) yields p(q, ui) (FT-LDOS). The loca- 
tion of the peaks of p(q, ui) provides information about 
the excitation spectrum of the electronic states^. In 
view of the suggestion of Ref. 2, it is particularly impor- 
tant to carefully examine how electron-boson coupling 
affects these peaks^. 

Modulations of the LDOS and the concomitant peaks 
in p(q, to) are due to electrons scattering from impuri- 
ties. In most studies the electron-impurity scattering has 
been treated either within a T-matrix approximation^ or 
a Born approximation 7 . In this framework, the peaks in 
p(q, u>) appear at wavevectors q = kf — k^, where kf and 
kj are points of high density of states at the energy ui. 
This approach has been quite successful in understanding 
FT-LDOS data for \u>\ < A^, and has been denoted as 
the 'octet' model^. We note that the success of this ap- 
proach depends on the existence of electronic states with 
long lifetimes in this energy range. Subsequent to this, 
there have been attempts to go beyond these approxima- 
tions by including the interaction of the electrons with 
either phonons or spin fluctuations^. 

The purpose of this paper is to reveal that for bias 



energies greater than ~ A^, where no quasiparticles are 
observed in photoemission^, there can be sharp peaks in 
p(q, to) due to a boson exchange process that appears be- 
yond the Born approximation for the electron-impurity 
scattering. For \ui\ > f2o; an electron can decay into a 
lower energy state by emitting a boson of energy Slo • As 
a result, the lifetime of the electrons with energy ui sig- 
nificantly larger than f^o is severely reduced, which is the 
primary reason for the absence of sharp peaks in p(q, uj) 
at the level of the Born approximation. But, going be- 
yond the Born approximation in the process shown in 
Fig. lb, an electron first emits a boson which reduces its 
energy to \oj\ — Slo . If this reduced energy is not signif- 
icantly larger than £7o, the resulting electron state (the 
internal fermion line of the diagram) is once again long- 
lived, and consequently contributes to peaks in p(q, uj) 
by scattering from the impurities. In this process, the lo- 
cation of the peaks at uj is determined by the excitation 
spectrum at an energy (\u\ — fio)sgn(w). That is, one 
has mirages of the octet peaks at \ujqi\ < A^ that are 
mirrored at an energy |o>qm| = \^>qi \ +^o- We note that 
this argument is based entirely on the energetics of the 
electron-boson interaction, and does not depend on the 
momentum space structure of the bosons (i.e., whether 
the bosons are spin fluctuations peaked at large wavevec- 
tors or phonons peaked at small wavevectors), though of 
course the momentum form factor of the bosons will lead 
to quantitative differences. 

Model. We study a two-dimensional system of super- 
conducting electrons described by a BCS model inter- 
acting with a boson mode, and coupled to an isotropic 
potential scatterer. It is described by the Hamilto- 
nian H = TLbcs + 'Him + Kel-b- Here Hbcs = 



,t 



-klCktJ 



where c 



(ckcr) creates (annihilates) electrons with spin a at 
wavevector k, the normal state dispersion is given by 
the tight binding expansion et = to + ^(cosfca; + 
cosfcj,)/2 + t2Cosk x cosk y + t 3 (cos2k x + cos2fc y )/2 + 
i4(cos 2k x cos k y + cos k x cos 2k y )/2 + t$ cos 2k x cos 2k y , 



2 



(a) 



k+q 



(b) x 
k p 



P+q k+q B 



(c) 




FIG. 1: (Color online) Panels (a)-(c) show the relevant dia- 
grams considered in this paper: Born diagram, vertex correc- 
tion, and the self-energy, respectively. The cross represents 
an impurity and the wiggly line a boson. In (a) and (b) the 
electrons (straight lines) are dressed by the self-energy. In (c) 
the electron line is bare. Panel (d) shows the constant energy 
contours of the bare dispersion for two energies: — lo — 20 
meV (the closed 'banana' contour) and — lo = 76 meV (the 
open contours), noting that A a = 37 meV. The wavevector 
qi connects the tips of the closed contour, and qA connects 
the inner of the two open contours along the antinodal direc- 
tion. 



and the superconducting gap has the d-wave form Ak = 
Ajvf(cosfc :r — cosk y )/2, with the lattice constant set to 
unity (M denoting the (ir,0) point). In order to study 
the sensitivity of p(q, lo) to the dispersion, we considered 
two sets of values for the parameters t{, one taken from 
Ref. [l^ and the other from Ref. Oj]. For both dispersions, 
the antinode is at = (l,0.18)7r. They differ in that 
the first dispersion has €m close to the Fermi energy (- 
34 meV), whereas it is further away for the second (-119 
meV). We focus here on results from the second disper- 
sion, though qualitatively similar results were obtained 
from the first. For Am, we choose 40 meV, a typical 
value for optimal doped Bi2212. 

The electron-impurity scattering is given by Tii m = 
X)k, q , CT c k+ qo - c k^ where Vq = leV in our calculation 
(note that Vq simply sets the scale for the FT-LDOS). For 
the sake of concreteness, we take the coupling between 
the bosons and the electrons to be of the form TC e i-b — 
9 £i ^ • where Si and Sj are the spin fluctuation and 
the electron spin operators respectively at site i, though 
our results hold equally well for optical phonons. We fix 
the magnitude of the coupling constant g from the condi- 
tion that in the normal state (Ak = 0) the inverse quasi- 
particle weight z _1 = 1 — 9R °^^ ~ 2 at the Fermi en- 
ergy, where is the electron self-energy due to inter- 
action with the spin fluctuations. This gives 3g 2 — 0.0176 
eV 2 . The dynamics of the spin fluctuations is given by 
X»u(Q,iftn) = x(*^n)<W withx(«fin) = 2^0/(^+^0), 
where x AIl /(Q,r) = (T t S, 1 (Q,t)S ii (-Q,0)) is the spin 
fluctuation propagator (the overall prefactor being ab- 
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FIG. 2: (Color online) (a) The real and imaginary parts of 
the 11 component of the self-energy, (b) Contribution to the 
FT-LDOS from the vertex diagram for — lo — Qo + Aa = 76 
meV. The dashed line shows the contribution from the Born 
diagram. The inset shows the contribution from the Born 
diagram for electrons without the self-energy correction. 



sorbed into the definition of g). Here /z, v are spatial 
indices, 0„ is a bosonic Matsubara frequency, and the 
mode energy f2n is taken to be 39 meV. In our model, 
we take the bosons to be independent of momentum for 
the following reasons. First, it allows us to concentrate 
on the energetics. As discussed in Ref. [HI as long as the 
form factor of the bosons is finite throughout the zone, 
then the self-energy will be dominated by the density 
of states singularities associated with the antinodal and 
M points of the internal fermion line. Second, it greatly 
simplifies the calculation of the vertex diagram (Fig. lb). 

Using Nambu notation, the 2x2 electron Green's func- 
tion with self-energy correction (Fig. lc) is given by 



G _1 (k,z) 



z - e k - £ii(z) -A k -£i 2 (k,z) 
-A k - EJ 2 (k,z) z + e k - £22(2:) 



where z is the complex frequency. The diagonal self- 
energies at T=0 are given byi£ £11/22 (*) = % 2 Skl z ^ 
e k {l + n /E k )]/[z 2 -{n + E k ) 2 }, where E k = + A^. 
In the above, the k-summation is performed numeri- 
cally with an intrinsic lifetime broadening factor r\ = 2 
meV for z = lo + irj. In Fig. 2a we show the real 
and imaginary components of Eu(z) as a function of 
lo. ImSn(cj) is zero for \lo\ < firj, & n d has pronounced 
peaks at \lo\ = Q + and ^0 + Em- Consequently, 
for energies up to and slightly beyond f2rj) the elec- 
tronic states form well defined quasiparticles, while for 
I a; I 3> r2o they do not. Strictly speaking, for momen- 
tum independent bosons, the off-diagonal self-energies 
vanish due to the ri-wave symmetry. However, in or- 
der to keep the antinode gap energy unrenormalized, 
we use the ansatz £12 (k, a;) = Ait[Z(w) — 1], where 
Z(lo) = 1 — ^jRe[En(w) + E 22 (w)]. In this scheme, the 
renormalized dispersion is E k m [A 2 . + e 2 ./ 'Z 2 (£; k )] 1 ^ 2 • 

Results. We first compute the contribution to the FT- 
LDOS from the Born diagram (Fig. la). This is given by 
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FIG. 3: (Color online) (a) Born contribution, and (b) Born 
plus vertex contributions, to the FT-LDOS for energies rang- 
ing from near — fio (-40 meV, top curve) to — Qq — Aa (-76 
meV, bottom curve) in steps of 2 meV. The curves are offset 
for clarity. 



pW(q,u>) = -(2/7r)Im[y E k Gf Q (k, W )(f 3 )^G^(k + 
q, u>)], where summation over repeated Nambu indices 
a, [3 — (1,2) is implied, and r» are Pauli matrices in 
Nambu space, with R denoting retarded propagators. 
The variation of p^(q, uj) versus q = (q x ,0) is shown 
as the dashed curve in Fig. 2b for — uj = A a + fi - We 
note that, due to large lifetime broadening, there is no 
peak in pW (q, uj) at this energy. The effect of the life- 
time can be seen clearly by comparing this curve with 
the one in the inset which is obtained by computing 
p( b \q.Lj = —fio — Aa) with £ = (i.e., for unrenor- 
malized electrons). 

Next we calculate the FT-LDOS contribution of the 
vertex diagram (Fig. lb). This is given by 



pM(q,L>) = - l-)Im 



^GfJk )W )T*( q ,a,) 



G*(k + q,a,)] 



(1) 



where 



T Q/3 (q,iw„) = — ^2 x{^n)G ai (p,iuj n - iCl n ) 
x (T 3 ) 7 sG S p{p + q,iuj n -in n ). 
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Here (3 is the inverse temperature and uj n is a fermionic 
Matsubara frequency. We note that in our model, the 



FIG. 4: (Color online) Dispersion of the the Fourier peaks 
along the bond direction versus energy for (a) the real part 
and (b) the imaginary part of the Born and vertex diagrams. 
Note the characteristic 'A' shape of the Born dispersion and 
its Qq (39 meV) displaced 'mirage' in the vertex dispersion. 



matrix T does not depend on k due to the momentum 
independence of the bosons. The variation of p( v ' (q, uj) 
along the bond direction is plotted in Fig. 2b for uj = 

— (f2 + Aa). We note the pronounced structure with 
a sharp maximum at q x — 0.307T and a sharp minimum 
at 0.34-7T. As \uj\ increases from uj = — (Slo + Aa) to 

— (f2o + Em), where Em is the renormalized energy at 
the M point (~ 60 meV), this structure evolves into a 
single broad minimum at q x — 0.237T (not shown). This 
can be contrasted with the Born result, where there is 
only weak structure in this entire energy range. 

For \ui\ ^> f2 , due to a large ImS aa (w), the k- 
summation in Eq. (1) involving the product of two G 
functions yields a quantity (denoted as C) which is ap- 
proximately constant as a function of q. The variation of 
p( v )(q, uj) with q for these energies is mainly due to that 
of T^p(q, uj). For the frequency summation involved in 
the computation of X^(q, uj), the main contribution is 
due to the bosonic pole which puts the electrons with 
momentum p and p + q (Fig. lb) at an energy uj + Qq 
for uj negative. This contribution can be written as (for 
uj < -n ) f R {q,uj) coh = 3g 2 J2 p G{p,UJ + n )T 3 G(p + 
q,uj + fio) - This expression is reminiscent of the Born 
contribution at an energy uj + f2o, and gives rise to sharp 
structure for \tu \ — £Iq < f2o since the electronic states are 
well defined at these energies. 

To understand this in greater detail, we plot in Fig. 3a 
the contribution to the FT-LDOS from the Born dia- 
gram, and in Fig. 3b the sum of the Born and vertex 
contributions, versus (q x ,0), for fio < \u>\ < Slo + Aa- 
In Fig. 4, we in turn plot the resulting peak dispersions 
from the real and imaginary parts of the Born and ver- 
tex diagrams for all |w| < f2o + Em- We note that the 
Born dispersion is well defined in the energy range be- 
tween and Em, although the Born peaks are damped 
for Oo < \uj\ < Em due to lifetime broadening. For 
|a->| < Aa, there are two peaks. The one at larger q x 
corresponds to scattering between the tips of the of the 
constant energy contours, which look like bananas in this 
energy range. This is denoted by the vector qi in Fig. Id 
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(structure corresponding to the even larger vector qs of 
Ref. is not plotted in Fig. 4). The one at smaller q x cor- 
responds to the so-called Tomasch peak noted in Ref. 0- 
It corresponds to where the banana first stops overlap- 
ping its q x displaced image. For < |cj| < Em, one 
finds a dominant maximum which traces out the separa- 
tion of the inner contours in Fig. Id along the antinodal 
((n, 0) — (ir, tt)) direction (denoted as in Fig. Id), with 
secondary peaks (not plotted in Fig. 4) corresponding to 
connecting an inner to an outer contour or an outer to 
an outer contour (these secondary peaks have less weight 
due to the reduced spectral weight of the outer contours). 
The combination of these peaks (two for \uj\ < A a and 
one for \u>\ > A a) gives a characteristic 'A' shape to 
the overall bond oriented dispersion, as is obvious from 
Fig. 4. Note there are some differences in the dispersions 
associated with the real and imaginary parts of the Born 
diagram and their connection to the vectors denoted in 
Fig. Id. This is due to several factors: the finite lifetime 
of the electronic states, the dispersion et, and the fact 
that ReG has a zero where ImG has a pole. 

We now turn to the vertex diagram. Its dispersion 
(Fig. 4) essentially mirrors the Born dispersion at a bias 
energy displaced by ^Iq. As a consequence, we denote this 
as a 'quasiparticle mirage'. In addition, for biases near 
ilo, the energy undisplaced Born term is also reflected 
in the vertex diagram (since the external lines in Fig. lb 
have well defined spectral peaks, and the boson exchange 
process has a limited phase space, for these energies). We 
note that there are some differences in the lineshapes of 
the imaginary part of the vertex term as compared to 
that of the Born term displaced by fio, as is evident in 
Fig. 3. This occurs since both the real and imaginary 
parts of the components of T contribute to p( v > (q, uo) as 
C is a complex quantity. 

Next, we comment on a few qualitative aspects of 
the vertex contribution to the FT-LDOS. (i) The in- 
clusion of a momentum form factor for the boson prop- 
agator should not make any qualitative change to the 
peak structure. Such a form factor (peaked around some 



Qo) can be thought of as a momentum constraint forc- 
ing k ss p + dQo (where a is a lattice group operation). 
However, in the mechanism discussed above, the elec- 
trons with momentum k and k + q (external lines of 
Fig. lb) do not play any special role, (ii) It is important 
that the boson that provides the 'step down' in energy 
has a sharp spectral function. A finite lifetime of the 
boson, or its dispersion with q, will broaden the Fourier 
peaks. For similar reasons, we anticipate that higher or- 
der vertex corrections will lead to weaker and broader 
contributions to the Fourier peaks because of the addi- 
tional momentum sums involved, (iii) Recently, the ob- 
servation of peaks in the Fourier transformed d 2 I/dV 2 
spectrum at q (0.4-7T, 0) has been reported for optimal 
doped Bi2212 at \w\ ~ Oo + Aa 2 ■ This is comparable to 
the peak position we find from our vertex corrected cal- 
culation at this bias energy. So far, though, no dispersion 
of these peaks has been reported. 

Conclusion. We demonstrated the importance of ver- 
tex corrections to electron impurity scattering in the 
study of tunneling spectroscopy data for the cuprate su- 
perconductors for absolute bias energies larger than f^o, 
where fio is the excitation energy of a boson coupled to 
the electrons. The vertex correction is due to emission 
and re-absorption of a boson by the electrons which leads 
to a 'step down' of the internal fermion line to energies 
where quasiparticle states are well defined, which as a 
consequence gives rise to sharp peaks in the FT-LDOS. 
We denote these new peaks as 'quasiparticle mirages', 
whose dispersion mirrors the previously observed quasi- 
particle interference peaks at absolute biases smaller than 
Hq. The observation of these dispersive 'mirages' would 
be an important reflection of the nature of the many- 
body interactions in cuprates. 
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